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Abstract 

The Kondo problem, for a quantum dot (QD), subjected to an 
external bias, is analyzed in the limit of infinite Coulomb repulsion by 
using a consistent equations of motion method based on a slave-boson 
Hamiltonian. Utilizing a strict perturbative solution in the leads-dot 
coupling, T, up to and orders, we calculate the QD spectral 
density and conductance, as well as the decoherent rate that drive 
the system from the strong to the weak coupling regime. Our results 
indicate that the weak coupling regime is reached for voltages larger 
than a few units of the Kondo temperature. 

PACS: 72.15.Qm, 73.23.Hk, 75.20.Hr 
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1 Introduction 



The transport of electrons through a quantum dot (QD) is a subject of in- 
creasing interest for its many apphcations in nanotechnology. Although there 
exists exact solutions for the ground state of the QD connected to metallic 
leads PP, no such methods are applicable to nonequilibrium systems. Typical 
cases are those of (a) a bias or (b) an electromagnetic radiation between the 
leads of a quantum well [2j. It has been suggested that, in case (a), two 
Kondo resonances appear, each one associated with every lead 0,111 • Some 
experimental evidence has been presented supporting this view j3] . Still, it is 
not well understood how inelastic (decoherent) processes modify the strong- 
coupling regime associated with the Kondo resonance. To our knowledge, 
there are only a few theoretical approaches to the problem. The non-crossing 
approximation (NCA) 0,101, being essentially numerical, can give analytical 
results only in the limit Ln{V/TK) » 1 (V being the applied bias and Tk 
the Kondo energy), the same limit in which renormalized perturbation the- 
ory (RPT) [7],jH] is valid. It is important, however, to have analytical results 
for the decoherence rate, not only to make theoretical predictions, but also 
to be able to interpret correctly the results of calculations done for more 
complicated situations, i.e., when time- dependent potentials are applied to 
the QD 0-^21 • A version of the nonequilibrium Kondo problem has been 
solved exactly to all orders in the applied bias V in 

In this article, we address the nonequilibrium Kondo problem for a QD 
in the limit of infinite Coulomb repulsion U, using an equation of motion 
(EOM) method based on a slave-boson Hamiltonian. By means of a strict 
perturbative solution in the leads-dot coupling, T, to orders and T^, we 
have been able to identify the decoherent processes that broaden the Kondo- 
level, as a function of the applied bias, and break the strong-coupling regime. 
The EOM method has been applied before to the Kondo problem |l4j but 
it was mainly used with approximations valid only for high temperatures 
, [ISj , [in] . We show that this method, when applied consistently to high 
orders, can be conveniently used in the nonequilibrium case, describing well 
the inelastic processes contributing to the Kondo resonance decoherence. We 
first check the method in the equilibrium regime, comparing our results with 
exact results using the numerical renormalization group (NRG) and also with 
the NCA. When used in the nonequilibrium case, we find that, for an applied 
bias larger than 2Tk, the system develops two-peaked resonant structures at 
±V^/2; we show, however, that decoherent processes take the system from 
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the strong-coupling regime, where the spin of the impurity is quenched by 
the lead electrons, to the weak-coupling regime, if V/Tk is larger than a few 
units. We conclude that our method provides an accurate and powerful tool 
to analyze nonequilibrium systems and time-dependent problems, situations 
where very few theoretical approaches are valid. 

2 Theory 

Our starting point for a QD level with infinite-U, coupled to two leads, is the 
slave-boson Hamiltonian for spin up and down electrons jl7j : 

H = Y. ^kaCl^ha + eo4^^'^ + Tkicl^Xa + xtcka) (l) 
ka cr PicF 

where the operators Xo- and are introduced in terms of an auxiliary boson 
field, b and b\ as Xo- = ^^d'o-, xj- = d\h and k = fc^, fc^j for left and right leads. 
In eq.((l}, the operator c].^ (cfco-) creates (anhilates) an electron of momentum 
k and spin a in one of the leads, while c/J. {d„) creates (anhilates) an electron 
of spin a in the dot. The constraint relation 

b^b + Ydida = l (2) 

(J 

is imposed to guarantee that no more than one electron occupies the QD 
since we are in the infinite-U limit. 

The physical (retarded) Green function for this problem is defined as 

mm 

« Mt); xiit') »= -I < t{t')xAt) + x^mlit') > d{t - 1') (3) 

and its Fourier transform is written as << Xa! >> ij^)- The EOM for this 
Green function reads (omitting the argument uj in all the Green functions): 

[w-eo) « Xa^xl »= 1- < > +J2'^k « (1 -n-a)cka;xt » + 

k 

YTk « d-adaCk-a,xi » (4) 

k 
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where the new Green functions appearing on the right-hand side of eq. (j3)) are 
defined as in eq.Q- n_fj is the number operator for electrons of spin —a in 
the QD and < A> will denote the expectation value of A from now on. In a 
second step, the EOM for these Green functions lead to a second generation 
of Green functions, viz., 

{w - ek) << Ck^]xl>>=Tk <<x„]xl>> (5) 



{W - tk) « d^^adaCk^a, xl »= " < X-a^k-a > + 
J^Tk' « Xcjcl'-a^k-a, Xl» + ^Tk' « X-aCk-aCk'a, xl » (6) 



{w - ek) « n_„Cka] xl »= J2 ^k' « X-aci_^Cka, xl » + 

k' 

J2 Tk' « X-aCk'^aCka, xt » (7) 
k' 

Note that the second generation of Green functions involve the excita- 
tion of two electrons in the leads. In our approximation to O(T^), we con- 
tract operators Cka, cl,^, where possible and obtain a closed equation for 
<< Xc7;xl » ■ In particular we take << Xacl^^^Ck-a] xl »= ^k,k' < 

flk-a >« Xa, xl » and << X-aCk-aCk'a, x\ » = « X-ctA' -a^ka] xl » = « 

X-aCk'-aCka',xl- >>= 0. Thls ylelds the following approximation to << 
Xa] x}j », which is exact to O(T^) (see Ref. [19J): 

«x.;xl»= ^ ^> (8) 

where 

So = E-^ (9) 

k ^ ~ 



We can obtain << Xa',Xa >> to 0{T ) by making all possible contrac- 
tions of operators Cka, c^k'a ^ third generation of Green functions. For 
instance, the EOM for << y^-aCk' -aCka] >> in eq.((Zj) is: 
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(u; + eo - efe - e^/) << X-^Ck'-aCka] xt »= Tk' « n^aCka] xt » - 

Tk' « dladaCk'-a;xi » « d^-adah' -ac\n ^Cka\ xl » " 

k" 

J^Tk" « (1 - na)cl„_^Ck'-aCka; Xi » (10) 
k" 

which is approximated by 

(a; + eo - Cfc - efc/) << X-aCk'-aCka, xl »= Tk' « n^^Cka, xl » - 
Tk' « d}_„d^Ck'-a] xl » +Tk < fik >« dl^^daCk'-^a] xl » - 

Tk' < hk' >« (1 - n„)cka; xl » (11) 

This procedure yields the Green functions already obtained in eqs.(I3]|7j) and 
a new one, << haCka] xl »i which has to be calculated as in eq.©. The 
first generation of Green functions are thus obtained as: 

{Uk - Sl(tl'fc)) << dl„daCk^a\ xl »= - < X~aCk-a > + 

rr ^ - -t / - ^Tk' «h„Ck'a\xl» /1o^ 
Tk < Uk >« Xa\ xl » -Tk <nk> 2^ z ■ (12) 

u, ^k — Cfc' + Co 



{Uk - S2(t^fc)) << n-aCka\ xl »= ^-ii^^k) « n^Cka, xl » + 

Tk' « dl^d^Ck'- 



Tk <hk>J2 



nXi » 



^k ~ ^k' + ^0 



(13) 



{Uk - T.2{uJk)) « n^Cka, xl »= - < Xlcka > +Tk < flk >« Xa] xl » + 

(~ \ - - -t / - ^ Tk' « n^^Ck'g] xl » n A\ 
S3(u;fc) << n_^Cfca; xl » -Tk < Uk > 2^ z ■ (14) 

^k — Cfc' + Co 

with ojk = uj — ek from now on. Eqs. ()12lll^ define a set of coupled inte- 
gral equations in the energy levels of the leads, e^, whose solution yields the 
Green functions to be introduced in eq.(@]). In these equations we have ne- 
glected expectation values of operators such as < xtx-o-cI-o-Cfccr >, which 
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give small contributions to << Xcr! >> to OiT"^) and have defined several 
self-energies given by: 

U^,) = T.T>A ~ ] , ^ + ~ .] ] (15) 

where vi = I- < >, 1^2 = I and 2/3 =< fik' >, (Sal^^fe) = ^2{'^k) - 
T,i{uJk)). These self-energies are associated with different Green functions 
and play a crutial role in the Kondo-resonance decoherence, as we will see 
below. Eqs. (|12llT^ reproduce the conventional T^-approximation of by 
neglecting all the terms proportional to T^. One can get, however, a new 
(and effective) T^-order near the Kondo resonance by neglecting the terms 
that couple these equations in an integral way. This approximation yields 
the following Green function: 



<< Xa;xt »-- 



[1- <n_,> - ^ — ^ . + 2^ Tk— 



t^fc-Si(u;fc) (u;^ - S2(c^fc))(u;fc - £4(0;^:))^ 

<nk> 

^ LJk-^iiu)k) ^ (cufc - S2(u;fc))(tl'fe - E4(a)fc))^ 



1^ - .0 - Eo - E "t,^^''^' J - (16) 



with 



£4(0;^) = T,2{uJk) + ^ ^^^'^'']^ (17) 

i^k — ^2{i^k) 

In eq.f|T?)|l. < > and < X-o-C/fc-o- > have to be calculated self-consistently 
At zero temperature and in a wide-band approximation, they are 

given by: 



-00 IT 

and 

r° c?^^^ Tfc « Xa;x l » (^ 

-00 vr LJk 
where 0; = is the Fermi level. 



duj + / N /X 

/ Im « Xa, xl » (18) 

J -00 71 
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3 Results 



In this section we first discuss some consequences that our approximation 
bears in the description of the Kondo peak. Notice that the T^-order of eq. (jH]) 
is given by eq. ()16|) neglecting the last terms of the numerator and denomi- 
nator and taking T,i{ujk) = 0. These new terms yield, however, contributions 
that are proportional to around the Kondo peak. This can be seen by 
considering the self-energies Si, S2 and S4, whose imaginary parts tend to 
broaden the Kondo-like peak associated with the factors J^k -^^-^^r^- In the 
T'^-order, eq. (fT5|) yields, for V = 0, eo < 0, |eo| >> A (2A being the one- 
electron level-width) and a wide-band approximation, S2 = 2So = —2iA, 
S3 ~ S2 and Si = 0, if \uj — ek\ << |eo|- Then, near the Fermi level, we 
obtain: 

D being the half-bandwidth, and 

The important point is that the term given by eq. (j2H) . which is neglected in 
the T^-order approximation, also gives a T^-contribution to << XajXl >> 
going like Ln{uj) at the Fermi energy. In particular, eq. lfT^ yields the Kondo- 
energy 

Tj, ^ D'/\AAy/'exp[^^] (22) 

intermediate between Dexp[—TT\eo\/ A] (the T^-order) and the exact result of 
Haldane {DAY/^exp[-Ti\€o\/2A] 

A check of our solution, eq. lfT^ . can be obtained by considering the case 
V = 0; this equilibrium case can be compared with the NCA and with 
the numerically exact calculations of Ref.|2I] using the NRG. Fig.l shows 
our calculated spectral density for the Kondo regime defined by eo = — 4A, 
A = 0.01 and D = 1 at zero temperature. Two cases, T^-order and eq. (|TT)|) 
are shown as a function of oj/Tk- Although for this particular case our value 
of Tx, eq. ()22j) . is very similar to Haldane's, we use this last one from now 
on. Our results indicate that eq. lfT^ yields a much better spectral density 
for the Kondo-resonance than the T^-order aproximation, with values very 



7 



similar to the ones given by the NRG above the Fermi energy, although 
our results yield a too small value of the spectral density below the Fermi 
energy; this deviation extends to energies of around 2Tk and suggests that 
a much more accurate solution is going to be obtained for voltages larger 
than this value (see below). Fig. 2 shows the same results for a mixed- valence 
regime (eg = 0, A = 0.01 and D = 1) aX zero temperature. In this case, 
eq. ()16|) yields an excellent approximation to the exact NRG-calculation and 
the Friedel-Langreth sum rule is satisfied with an accuracy better than 10%. 
The agreement is even better in the empty orbital regime. Fig. 3, where the 
sum rule is fulfilled to within 5% accuracy. We should mention that our 
present results are in much better agreement with the exact results than the 
ones obtained using Lacroix's approximation J3] and shown in [T^ . 

Now we turn our attention to the nonequilibrium case V ^ Q 'm. the sym- 
metric situation, where the chemical potentials of the left and right leads are 
±1^/2. A fully consistent solution of the case V ^ Q requires the application 
of the Keldysh formalism, as explained in ^H] ■ The symmetric case can be 
analyzed, however, at the level of the above equations because, to good accu- 
racy, < > and < X-o-Cfc_o- > do not vary with V for the values used here. 
We start by analyzing the self-energies Ei and S4 in the regime ujk — >■ 0, 
which is the relevant regime for analyzing the Kondo-resonance decoherence. 
It is easy to see that eq. (fTH|) produces — » 0) , as in the equilibrium 
case and therefore, we find no decoherence for V 7^ 0. It is interesting to 
proceed further with the EOM method and see how these self-energies are 
modified. By going to O(T^), we find for Ei and S2: 





^'^cDfe-efc/+e_(/c')+iA Cjk + - e+{k') + il^ 
T2 (1- < hk' >) 1 Tl„{l- < Uk" >) 



+ 



t^fc - Cfc' + e- (^0 +iAuJk + ek" - tOk + ^k" - e+ (^") + 
Tf,(l- < hk' >) 1 Tl„{l- < Uk" >) 



X 




+ 0(r*) (23) 



i^k + Cfc" - eg t^fc - eg + e_(g) + iA 



and 



E,{^k) = T.T'k'[ 



1 1 



] + 



k' 



^k — e/c' + ^-{k') + iA uJk + ek' — e+(/c') + iA 
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^„ oJk - ^k' + (^-{k') + iA LUk + tk"-eki Uk + ^k" - (^+{k") + 

where e±(efc/) = cq - f Ln(^^^^^). We also get T,3{uJk) = ^2{<^k) - ^li^^k), 
T,Jujk) = T.2(tOk) + -^J^l^'}- N , as before. These higher order self-energies give 
rise to the first non-zero contribution to the Kondo resonance decoherence, 
due to the cross terms between the right and left leads, as has already been 
pointed out in [3]. Eqs. ()23|) and (j^^ represents the first terms of an expansion 
for El or S2, using (T^/eo)^ as the expansion parameter, and allows us to 
analyze the Kondo resonance decoherence as a function of an external bias. 
For = 0, we again find T,i{uk ^ 0) ^ and S4(tl'fc ^ 0) ^ but for 
V ^ these self-energies develop no n- negligible imaginary parts for cu^ — * 0. 
In the bias regime, defined by V << |eo| < A, we find 

0; V) ~ - 2^Ln^) = -27 (25) 

Co 2tt TT\eo\ |eo| 

and S2 — —22 A. In eq. ()25|) . eo is a renormalized level given approximately 
by eo — Co + ^Ln^ . Fig.4 shows the evolution of the Kondo-like resonance 
spectral density for the same parameters used in Fig.l, but introducing dif- 
ferent biases on the QD. Here we have used eq.()16p and the self-energies 
to T^-order. Our results are compared with the NCA-calculations of Ref. 
jHj, where available. The agreement between both approaches, for > 0, 
is remarkable. The difference between both cases can be attributed to the 
bandwidth used in Pj, which is ten times smaller than ours. For a; < 0, how- 
ever, the agreement is only reasonable at low bias (like in the equilibrium 
case shown in Fig.l), which improves as V increases. 

We have also calculated the QD-conductance which, in the symmetric 
configuration we are considering, is basically proportional to the mean value 
of the spectral density in the interval (—1^/2, V/2). Fig. 5 shows the conduc- 
tance G in units of the conductance quantum Go = 2e^//i, as a function of 
V/Tk] a good fit is obtained by the expression: 

G/Go = 1/(1 + 0.75Ln(l + {^f)), (26) 

for large values of V/Tk- These values also show good agreement with the 
NCA-calculations of [7.^ for V > 2Tk, where our solution can be taken with 
great confidence, although in this case, G/Gq behaves like l/Ln^CV/TK) for 
very high voltages, Ln(y/TK) » 1. 
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In addition, we have also addressed the important issue of strong versus 
weak couphng regimes, as a function of the bias. It has been argued 
that the strong-couphng regime might subsist, even for high bias {V » 
Tk), if the decoherence rate 7 is smaller than T* = —V/2 + \Jt'^ + {V/2Y. 
With 7 given by eq.(j2Sl), the strong-coupling regime is obtained for V < 
vr^/^^Tft-. This is a very stringent condition; for instance, for the values 
of the parameters in Figs. la and 2, we find strong coupling for V < 3.5Tk, 
for which the conductance is larger than 0.4Go. We conclude that, when 
subjected to an external bias V > T^, the QD-Kondo resonance develops a 
two-peaked structure at ±V/2 for V > 2Tk; the system, however, evolves 
quickly as a function of the voltage, moving from a strong-coupling regime 
with the spin impurity screened by the electrons of the leads, to a weak- 
coupling regime with that spin unscreened. 

4 Conclusions 

The EOM method provides an accurate approximation to the Kondo problem 
in and out of equilibrium situations, the only requirement being to include 
all the Green functions in a strict expansion in T^. Our results to T^-order 
in equilibrium provide a good description of the Kondo peak in the Kondo 
regime and an excellent aproximation to (numerically) exact results in the 
mixed-valence and empty-orbital regimes. In the nonequilibrium regime, it 
is necessary to carry the method to O(T^) to find the Kondo resonance deco- 
herence, since it involves the coupling between the right and left electrodes 
through the QD. In particular, in our approach, the decoherence rate is found 
analytically to be nearly proportional to V, which indicates that the system 
is driven readily to the weak couping regime ifV> 7rV2MT^, with the spm 
impurity unquenched by the electrons of the leads. 
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Figure 1: The spectral density in equilibrium as a function of w/Tk for 
D = 1 and A = 0.01 is plotted for the Kondo regime: eo = — 4A, T^/A = 
1.87 X 10-2. Dotted line: T^-EOM; full line: T^-EOM; dashed line: NRG 
from [21]: dashed-dotted line: NCA from [21] (T = IQ-^D). 
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Figure 2: The spectral density in equilibrium as a function of w/Tk for 
D = 1 and A = 0.01 is plotted for the mixed- valence regime eg = 0, T^/A = 
1. Dotted line: T^-EOM; full line: T^-EOM; dashed line: NRG from [21]: 
dashed-dotted line: NCA from [21] (T = IQ"^D). 
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Figure 3: The spectral density in equilibrium as a function of w/Tk for 
D = 1 and A = 0.01 is plotted for the empty-orbital regime eo = 2A, 
Tk/^q = 1. Dotted line: T^-EOM; full line: T^-EOM; dashed line: NRG 
from [21]: dashed-dotted line: NCA from [21] (T = IQ-^D). 
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Figure 4: The spectral density in the nonequihbrium regime as a function 
of w/Tki calculated with our EOM method for the same parameters as in 
Fig.l is shown for: a) V = 2Tk (continuous line) and V = 6Tk (dotted line). 
In Figs. 4b, 4c and 4d our results (continuous lines) are compared with the 
NCA results of Ref.[10] (dashed lines) for: h) V = 8Tk, c) V = 20Tk, d) 
V = AOTk. 
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Figure 5: The conductance (in units of the quantum of conductance) calcu- 
lated for the same parameters as in Figs.l and 4 (continuous line) is compared 
with the fitting function of eq.(|^ (dashed line). 
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